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ABSTRACT 

We analyze the confining effect of cloud collapse on an expanding supernova shockfront. We solve 
the differential equation for the forces on the shockfront due to ram pressure, supernova energy, and 
gravity. We find that the expansion of the shockfront is slowed and in fact reversed by the collapsing 
cloud. Including radiative losses and a potential time lag between supernova explosion and cloud 
collapse shows that the expansion is reversed at smaller distances as compared to the non-radiative 
case. We also consider the case of multiple supernova explosions at the center of a collapsing cloud. For 
instance, if we scale our self-similar solution to a single supernova of energy 10 51 ergs occurring when 
a cloud of initial density 10 2 H/cm 3 has collapsed by 50%, we find that the shockfront is confined to 
~ 15 pc in ~ 1 Myrs. Our calculations are pertinent to the observed unusually compact non-thermal 
radio emission in blue compact dwarf galaxies (BCDs). More generally, we demonstrate the potential 
of a collapsing cloud to confine supernovae, thereby explaining how dwarf galaxies would exist beyond 
their first generation of star formation. 

Subject headings: galaxies: individual (SBS 0335-052, Henize 2-10) - galaxies: dwarf - supernova 
remnant - galaxies: starburst - stars: dwarf nova 



1. INTRODUCTION 

Blue Compact Dwarf galaxies (BCDs) are observed to 
be experiencing intense star formation in a spatially com- 
pact region. The timescale for this star formation episode 
of mass ~ 1O 8 M is only a few million years. Massive 
stars have short lifetimes and thus supernova (SN) ex- 
plosions may occur concurrently in regions of cloud col- 
lapse. It is quite probable that this occurs in blue com- 
pact dwarf galaxies. In this paper, we study the confin- 
ing effect of cloud collapse on an expanding supernova 
shockfront. 

The forces acting on the shockfront include those due 
to supernova energy, ram pressure and gravity. When 
a supernova shockfront expands into a stationary cloud 
(Sedov 1946, Taylor 1950, Bisnovatyi-Kogan & Silich 
1995), the accreted mass does not contribute to a change 
in the momentum flux. In a region of cloud collapse, not 
only is mass accreted at a faster rate, it also contributes 
to a change in the momentum flux. This additional ram 
pressure leads to the confinement of a supernova shock- 
front in a collapsing cloud. The potential of ram pres- 
sure to confine stellar winds of a noncentral OB star to 
generate a steady-flow situation has been shown (Dopita 
1981). We demonstrate that ram pressure can also con- 
fine expanding supernova shockfronts with pressures as 
high as ~ 10~ 9 dynes/cm 2 as compared to stellar wind 
pressure which is ~ 10~ 12 dynes/cm 2 (Dopita et al 1981). 
Our self-similar solution describes the time evolution of 
a central supernova explosion in a collapsing cloud. 

We start with a discussion of the gravitational collapse 
of a cloud, and then obtain the equation for the expan- 
sion of a SN shock in the cloud. We first neglect radiative 
losses and then include them. Next, we consider a time 
lag between when the supernova explodes and when the 
cloud begins to collapse. Finally, we consider the case 
of multiple SN explosions occurring at the center of a 
two-component collapsing cloud. The inner component 
has higher density and collapses rapidly to produce the 



central stellar population. A few million years later, the 
resulting supernova shockfronts collide with the infalling, 
less dense, outer component of the cloud. 

2. MODEL 

We first consider a simple model for the free-fall gravi- 
tational collapse of a cloud with negligible pressure. The 
gravitational force on a differential spherical shell with 
initial radius i?o and radius R(t), is 

SR _ _ GM(Rg)dM{R) 

dt z R z 
where G is the universal gravitational constant and 
M(Rq) is the initial mass enclosed by this shell. 
M(Rq) = AttRqPo/3, where po is the initial density of 
the cloud which is assumed to be uniform. Solving this 
equation, assuming the initial velocity is zero, gives 
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The free-fall time of the cloud is 
_ f R ° dR ( 
tff ~J -dR/dt ~ \32Gp Q 

«5.2 x 10V (^r 172 , 

where is density in units of 10 2 H/cm 3 . 
fall time is independent of the initial radius, 
collapse is homologous as shown in Figure 1. 

The velocity of the shell hitting the shock front con- 
tributes to the ram pressure that slows down the front. 
We normalize the time and radius as r = t/tff and 
R = R/Ro. Hence we obtain 

This integral gives a cumbersome expression for R(t). 
Therefore, we approximate the dependence as 
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Thus the velocity of the cloud at the shockfront radius, 
Rsh, is 

wRshil-fir)} 1 / 2 
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Accretion of mass just outside the shock front acts to 
increase in density. Mass conservation gives 

AM = AnRlARopo = 4nR 2 ARp(t) . (7) 

Using R = R f{r) and AR = Ai? /(r) gives 
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There are three forces experienced by the shock front: 

(i) The energy released by the supernova explosion Esn 
inside the shell gives an outward directed pressure p = 
(7 - l)E SN /(±irR 3 sh ) or force ^R 2 sh p. We use 7 = 5/3. 

(ii) The ram pressure of the cloud shell hitting the shock 
front results in an inward directed force (or pressure). 
When the direction of motion of the external medium and 
the shockfront are the same, the ram pressure acts only if 
the speed of the medium is greater than that of the shock 
front. This is accounted for by a Heaviside function, 
H , of the difference between the shockfront velocity and 
cloud velocity (with H(> 0) = 1 and H{< 0) = 0). (iii) 
The gravitational force is of course inward. 

The ambient pressure of the interstellar medium is ~ 
10 -10 dynes/cm 2 for a density ~ 100 H/cm 3 and tem- 
perature ~ 10 4 K. For a typical speed ~ 100 km/s, the 
ram pressure is ~ 10 -8 dynes/cm 2 . For this reason, we 
neglect the ambient pressure of the interstellar medium. 

The rate at which the SN shockfront accretes mass is 
proportional to the relative velocity between the cloud 
and the shockfront. Thus, 



= AirR sh p U 



where 
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The accreted mass initially moved at the cloud velocity. 
It thus contributes to a net change in the momentum flux. 
For the case where the SN explosion and cloud collapse 
begin at the same time we have 
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We normalize time by the free-fall time and the shock- 
front radius by the Sedov- Taylor radius at the free- 
fall time (within a factor of 1.24); that is, R s h = 
R s h/RsT(tff), where 

RsT(t ff ) = (h-l)EsNt 2 ff /po) 1/5 ^52pc {E^) 1 ' 5 ^)- 2 ' 5 

_ (12) 
We normalize mass by M s h — M s h/[4:TrpoR ST (tff)/3]. 
We define a 'momentum' variable 
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Rewriting equation J!|J| and equation (|11|) in terms of 
the normalized parameters gives 
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Thus, we have a system of three non-linear coupled first 
order equations. For the initial conditions we assume at 
an early time, R s h and dR s h/dr are given by the stan- 
dard Sedov- Taylor solution for a stationary cloud. We 
use the initial condition that at r = 10 -6 , R s h — 0.0049, 
M sh = 1.2 x 10~ 7 and V= 0.00024. Figure 2 shows 
R s h(r) and the corresponding Sedov- Taylor solution. 
The infalling cloud at first slows the shock expansion 
and later reverses its motion. 

3. INFLUENCE OF RADIATIVE LOSSES 

In the above analysis we assumed that the energy 
within the SN shock was constant, but this is valid only 
for relatively short times. We next take into account 
the radiative losses using the pressure driven snowplow 
model of Cioffi, McKee, and Bertschinger (1988). We 
find that the supernova is confined to a much smaller 
region in a shorter time. We use the mean pressure p rl 
including radiative losses derived by Cioffi et al. (1988) 
which can be written as 
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where magnitude of a is 2.3 x 10 -17 and tjm is metallicity 
in units of solar metallicity. Thus p ri replaces our earlier 
pressure p = Esn / '{^R^h)- Next, we assume that the 
metallicity is solar metallicity and we use equation JHJ for 
density as a function of time. The differential equation 
including the radiative losses is 

4 ( M Kh ^\ - 4n R 2 sh (p rl - P H(U) U V cl )- GM ° h 
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We continue to normalize time by the free-fall time. How- 
ever, we now normalize the radius by 
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We again obtain a self-similar differential equation and 
we solve it using the initial condition that the shockfront 
follows the Sedov solution at r = 10~ 6 . We use this ini- 
tial condition because the transition from the Sedov solu- 
tion to the pressure driven snowplow radiative shockfront 
occurs at 



tpDS ~ 1-33 x 10 4 yr 
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(Cioffi et al. 1988). For E = 10 51 ergs and n = 
10 2 H/cm 3 , t PDS ~ 1700 yr ~ i///3000. 

The SN explosion may occur at a time At after the 
cloud collapse has begun. The velocity and density of 
the shells hitting the shock front is then larger. We take 
the zero of time t or r = t/t// to be the time when the 
supernova explodes so that only equation ® needs to be 
modified and replaced by /(r) = R = [1 - (r + At) 2 ] 1 / 2 , 
where At — At/tff is termed the "lag." 

Figure 3 shows sample results for different lags which 
shows that the shock expansion is more strongly deceler- 
ated the longer the lag. For no time lag, the shockfront 
is confined to B, s h = 1-99 and r = 0.51. In comparison, 
for 50% time lag, the shockfront is confined to R s h = 
1.33 and r = 0.23 and for 80% time lag, the shockfront 
is confined to R s h = 0.75 and r = 0.09. 

4. MULTIPLE SUPERNOVA EXPLOSIONS 

We now consider a more complete physical picture. 
Suppose that there is a localized clump of density 10 4 
H/cm 3 in a cloud of density 10 2 H/cm 3 . Suppose that 
the ambient outer cloud and the inner clump begin to 
free-fall. The clump collapses in ~ 0.5 Myrs and forms 
supcrnovae in another 2 — 3.5 Myrs (lifetime of 30 — 24M Q 
star). The free-fall time of the cloud is ~ 5 Myrs and 
thus, the time lag between when the supernova explodes 
and when the cloud begins to collapse ranges from O.Mff 
toOMff. 

When the localized clump collapses, the initial mass 
function is such that many low mass (lm) stars form 
at the same time as a few high mass stars. The dis- 
tribution is dN/dM oc M~ 2 - 35 (Salpeter 1955). If we 
take the high mass (hm) threshold at > 20M Q , then 
Mim ~ 2OOiV/ im M0. The mass of the low mass stars 
contributes to the gravitational force in addition to the 
forces described in equation 1)1 8(1 . The right-hand side of 
this equation becomes 
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We find however that the gravitational force contribu- 
tion of the low mass stars does not significantly alter the 
trajectory of the SN shockfront. 

5. DISCUSSION 

Regarding the stability of the shockfront, notice that 
the Rayleigh- Taylor instability occurs in a stratified re- 
gion where the effective gravity points in the direction of 
the less dense medium. In the case of an expanding su- 
pernova shockfront, the effective gravity in the frame of 
the shockfront is the shockfront acceleration. Our results 
for the trajectory are all concave down implying that ef- 
fective gravity points radially outward and thus, away 



from the relatively rarer medium. Hence, the shockfront 
is Rayleigh- Taylor stable. 

We conclude that ram pressure from a collapsing cloud 
reduces the shockfront velocity of one radiative super- 
nova explosion (occurring in a localized region when the 
ambient cloud has collapsed by 50%) to zero at a normal- 
ized radius of 1.33 and a normalized time of 0.23. These 
fractions scale by the energy of the supernova and the 
density of the interstellar medium. 

Observations of nonthermal radio emission in galaxies 
suggest whether or not the supernovae are confined. For 
example, a 22 GHz map of M82 extending over a few 
hundred parsecs is estimated to have 40 unconfined su- 
pernovae (Golla et al 1996). An 18 cm nonthermal radio 
emission map of VIIZwl9 extending over 310 pc is es- 
timated to have 2500 unconfined supernovae remnants 
(Beck, et al 2004). Two examples of apparently con- 
fined emission, SBS-0335 and Henize 2-10, are described 
below. 

In a model of the nonthermal radio emission from SBS- 
0335 the emission is confined to 17 pc (Hunt et al. 2004) 
within a 520 pc region (Thuan et al. 1997) where star 
formation occurs in six super-star clusters with ages < 25 
Myrs. If a supernova explosion of energy 10 51 ergs occurs 
after a cloud of initial density 10 2 H/cm 3 has collapsed 
by 50%, it is confined to ~ 15 pc in ~ 1 Myr. If this were 
a supernova exploding in a stationary cloud, it would ex- 
pand beyond ~ 20 pc in as little as 0.4 Myrs. This would 
destroy the cloud and the case of a young supernova rem- 
nant in a stationary cloud is considerably less probable 
than a supernova confined in a collapsing cloud. 

Very Large Array imaging of another BCD, Henize 2- 
10, indicates a < 8 pc region of 1 mJy radio sources in the 
central 5 starburst region. Henize 2-10 has HII regions 
of sizes between 3 pc and 8 pc and densities between 
1500 and 5000 H/cm 3 (Kobulnicky, Johnson 1999). If 
we hypothesize that a cloud of an average initial density 
of 3000 H/cm 3 has collapsed by 50% when 10 supernovae 
explode, we find that the supernovae are confined to ~ 7 
pc in ~ 0.2 Myr. Thus, cloud collapse successfully con- 
fines supernova explosions and can account for observed 
compact nonthermal radio emission. This simple model 
can help understand continuous or second/third genera- 
tion star formation since it suggests why the cloud is not 
devastated by first generation supernovae. 
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Fig. 1. — The collapse of cloud shells with different initial radii are shown as solid lines. The radius is normalized by the initial radius of 
the outermost cloud shell and the time is normalized by the free-fall time. The dotted line indicates the approximation of the trajectory 
used here. The approximation has a maximum deviation of 5.3% and a mean deviation of 2.6%. 
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Fig. 2. — Shock wave expansion R a h(T~) in a collapsing cloud and the standard Scdov- Taylor expansion in a stationary medium. Time 
T is normalized by the free-fall time and the radius is normalized by the Sedov- Taylor radius at the free-fall time. The shock expansion is 
slowed down and subcquently turned around by the collapsing cloud. 
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Fig. 3. — Shockfront radii including influence of the radiative losses. We also show the shockfront radius for a time lag of 50% and 80% 
of the free-fall time between supernova explosion and cloud collapse. Time zero is taken as the instant the supernova exploded. 
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Fig. 4. — Evolution of the shockfront radius for single and multiple SN explosions. 



